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Abstract 

GNU TEXmacs is a free Wysiwyg word processor providing an excellent 
typesetting quality of texts and formulae. It can also be used as an interface 



cn 

U 

ZJ ', to Computer Algebra Systems (CASs). In the present work, interfaces to 

three general-purpose CASs have been implemented. 

(N 
> 

^ \ 1 TfeXmacs 

t^^ i GNU T[T;Xmacs []I| is a free (GPL) word processor which 

• typesets texts and mathematical formulae with very high quality (like BTgX), 
Q ! • emphasizes the logical structure of a document rather than its appearance (like 

r\ \ • is easy to use and intuitive (like typical Wysiwyg word processors), 

• can be extended by a powerful programming language (like Emacs), 

• can include PostScript figures (as well as other figures which can be converted 
to PostScript), 

• can export ETgX, and import BTgX and html, 

• supports a number of languages based on Latin and Cyrillic alphabets. 

It uses T[t;X fonts both on screen and when printing documents. Therefore, it is 
truly Wysiwyg, with equally good quality of on-screen and printed documents (in 
contrast to LyX which uses X fonts on screen and calls I^T[t;X for printing). There 
is a similar commercial program called Scientific Workplace (for Windows). 

TEXmacs can also be used as an interface to any CAS which can generate I^T^X 
output. It renders I^T[t;X formulae on the fly, producing CAS output with highest 



typesetting quality (better than, e.g., Mathematica, which uses fixed- width fonts 
for formula output). A user can utilize editing facilities of TEXmacs: copy (a part 
of) a previous input into the new one, edit it and sent to the CAS, copy a result 
derived using the CAS into a paper, etc. In the present talk, I give some examples 
of using Maxima, MuPAD and REDUCE via T[T;Xmacs. It is not my aim to describe 
these powerful and complex CASs; I only show examples of typesetting produced 
by TgXmacs. This talk has been written in TEKmacs and exported to I^TgX. 

2 Maxima 

Macsyma is one of the oldest and most mature CASs. It was developed at MIT 
during the end of sixties - beginning of seventies. Later, it was owned by various 
commercial companies. Now it seems practically dead. 

Fortunately, a free CAS Maxima is now under GPL. It is based upon the Macsyma 
code base from seventies, with a number of later enhancements. It incorporates a 
lot of mathematical knowledge, is stable and well tested. From its very beginning, 
Macsyma (and Maxima) pays much attention to mathematical correctness; for ex- 
ample, if the form of an integral depends on the sign of a parameter, it will ask 
the user about it, or use an assumption - other systems only recently incorporated 
similar facilities. It is an excellent platform for research projects, because it provides 
a solid foundation, and it cannot vanish into thin air as commercial systems (e.g., 
Macsyma) can do at any moment. Its text-based interface now looks somewhat old- 
fashioned. When combined with a nice graphical interface provided by TgXmacs, it 
can compete with commercial CASs like Matliematica and Maple, and even produce 
higher-quality output. And it is difficult for commercial vendors to beat the price 

Here is a sample Maxima session within T^Kmacs. 

GCL (GNU Common Lisp) Version(2.4.0) Tue May 15 15:03:11 NOVST 2001 

Licensed under GNU Library General Public License 

Contains Enhancements by W. Schelter 

Maxima 5.6 Tue May 15 15:03:08 NOVST 2001 

(with enhancements by W. Schelter) . 

Licensed under the GNU Public License (see file COPYING) 

(CI) (x"2-y"2) / (x~2+y~2) +sin(alpha) "2 ; 

(Dl) ^~^. + sin^ a 

(C2) expand((x+y-l)"5) ; 

(D2) y'' + 5x?/^ - by^ + IQx^y^ - 20xy^ + lOy^ + lOx^y^ - SOx'^y^ + 30xy^ - lOy^ 

+ bx^y - 20x^y + 30x'^y - 20xy + 5y + x^ - 5x^ + Wx^ - lOx^ -|- 5a; - 1 



(C3) solve (a* x~2+b*x+ c , x) ; 

Vb^ -4ac + b 
(D3; x = :: ,x 



b-Vb^- ^ac 



2a 2a 

(C4) integrate (sqrt (x~2+a) ,x) ; 

Is a positive or negative? negative; 

_ ^ a\og (2\/x^~+~a + 2x) x\Jx^ + a 

(D4) ^^ H 

2 2 

(C5) assume (a>0); 

(D5) [a > 0] 

(C6) integrate (sqrt (x~2+a) ,x) ; 

aASINH (^) ^ 



\Jx'^ + a 



(D6) 

2 2 

(C7) integrate(exp(sin(x)) ,x,0,7opi) ; 



(D7) 



dx 



(C8) diff (f (x),x,2); 

(KB) |j/W 

(C9) g : gamma ( 1+x) ; 

(D9) r(x + l) 

(CIO) taylor(g,x,0,3); 

(672 + 7r2)x2 (273 + 7r27 + 4C(3))x3 
7X+ - 



(DIO) 1 



12 



(Cll) m:entermatrix(2,2) ; 

Is the matrix 1. Diagonal 2. 

4. General 

Answer 1, 2, 3 or 4 : 4; 



12 



Symmetric 3. Antisymmetric 



Row 1 Column 1 


a 




Row 1 Column 2 


b 




Row 2 Column 1 


c 




Row 2 Column 2 


d 




Matrix entered 




\c d ) 




(C12) m-~(-l); 
/ d 


h 


(D12) "'^-l" 


ad— be 
a 



\ ad— be ad— be 

(C13) fac(n):=if n=0 then 1 else n*fac(n-l); 

(D13) fac (n) := if n = then 1 else niac {n — 1) 

(C14) trace(fac); 

(D14) [fac] 



(C15) fac(5); 
1 Enter fac [5] 
2 Enter fac [4] 
3 Enter fac [3] 
4 Enter fac [2] 
5 Enter fac [1] 
6 Enter fac [0] 
6 Exit fac 1 
5 Exit fac 1 
4 Exit fac 2 
3 Exit fac 6 
2 Exit fac 24 
1 Exit fac 120 
(D15) 120 

(C16) f (x) :=block([a,z:0] ,a:x+l ,a:a/z,a+l) ; 
(D16) / (x) := block ([a, z : 0] ,a : x + l,a : -,a + l) 

(C17) debugmode(true) ; 

(D17) true 

(CIS) f(u); 

Division by 

— an error. Entering the Maxima Debugger dbm f (x=u) 

(dbm:l) a; 

u + 1 

(dbm:l) z; 



(dbm:l) :q 

(C19) plot2d(sin(x)/x, [x,-10,10]); 

(D19) 

(C20) f (x,y) :=sin(sqrt(x~2+y~2))/sqrt(x"2+y~2) ; 

(C21) plot3d(f (x,y), [x,-10,10] , [y,-10,10]) ; 
(D21) 

Plots appear in separate windows (Fig. |I]). When the mouse is over such a 
window, its coordinates are continuously displayed at the upper left corner (in the 
3d case, also z of the surface at the mouse position {x,y) is shown). When the 
mouse is at the upper left corner, a menu appears. It allows the user to control the 
plot: zoom, rotate (in the 3d case), print, save as PostScript, etc. 

The toolbar icon showing the question mark shows the Maxima documentation. 
The Maxima manual is in html; it is imported into TgXmacs and shown in a new 
buffer. Hyperlinks work with double click. It is easy to arrange things in such a way 
that you do calculations in one T[7iXmacs window, and read the manual in another 



Plot of z ■■ 





Figure 1: Maxima plots 
window. 

3 MuPAD 

MuPAD [Q is the most recent addition to the family of universal CASs. It is being 
developed at University of Paderborn, and commercially distributed by SciFace. In 
some cases, it can be obtained free of charge (see their web site), but it is not free 
software. It is rather similar to Maple, but designed from scratch, and some new 
fundamental ideas were incorporated. The library is not so extensive as those of 
older systems, but MuPAD is progressing fast. Its interface is text-based (except 
the Windows version). Therefore, adding a high-quality graphical formula output 
provided by TgXmacs is useful. 

Here is a sample MuPAD session within TgXmacs. 

* * MuPAD 2.0.0 — The Open Computer Algebra System 

/I /I 

* * I Copyright (c) 1997 - 2000 by SciFace Software 

I * — I -* All rights reserved. 

1/ 1/ 

* * Licensed to: Andrey Grozin 



> (x"2-y"2)/(x~2+y~2)+sin(alpha)~2 

sm a +^^- — ^ 



:^ expand ((x+y- 


1)~5) 


















5x + 5y - 


■20xy 


-10a;2 


+ 10x3- 


-lOy' 


-5x^ 


+ 10y' 


+ x' 


-5y' + y' 


+ 30 


xy 


+ 30x'^y 


-20xy 


3-20 


x^y + 5x 


y' + 5 


4 

X y- 


SOx'y' 


+ 10 


x^ y^ + 10 


x'y'- 


- 1 



> solve (a*x~2+b*x+c=0,x) 



C ifa = OAb = OAc = 

{} ifa = 0A6 = 0Ac^0 

{-f} ifa = 0A6^0 



h Vb-^— 4ac y/b^ —i a c b 

'2 2 2 2 

a ' a 

> int(sqrt(x"2+a) ,x) 

X a/oTo? a In (z + Va + x"^) 



if a^O 



2 2 

> il:=int(exp(sin(x)) ,x=0. .PI) ; float(il) 

r 

exp (sin (x)) dx 

6.208758036 

> diff (f (x),x,x) 

> g:=gainma(l+x) 
7(2; + l) 
^ series(g,x=0,4) 

' Vl2 2 y V 3 6 12 ' ^ ^ 

> M:=matrix([[a,b] , [c,d]]); 1/M 
a b 
c d 



-ad hc~ad 



he— ad bc~ad 

> plotfunc2d(sin(x)/x,x=-10. .10) 

Warning: Dumb terminal: Plot data saved in binary file save.mp 
[plot] ; during evaluation of 'plot2d' 

> f :=(x,y)->sin(sqrt(x"2+y"2))/sqrt(x"2+y~2) 
(x, y) ->sin(sqrt(x"2+ y"2))/sqrt(x~2+ y"2) 

> plotfunc3d(f (x,y) ,x=-10. .10,y=-10. .10) 

Warning: Dumb terminal: Plot data saved in binary file save.mp 
[plot] ; during evaluation of 'plotSd' 

> quit 
The end 

Plots appear in separate windows (Fig. ^). Spurious warnings about dumb ter- 
minal may be ignored. Plots are displayed by the program vcam which is distributed 
with MuPAD. They can be controlled (and saved to PostScript files) via menus. 

The question mark icon opens the help menu. MuPAD documentation is exten- 
sive and high-quality. It is written in an extension of I^TgX with hyperlinks, and 





Figure 2: MuPAD plots 

displayed by an extended dvi viewer distributed with MuPAD. Therefore, it has 
high typesetting quality. 

4 REDUCE 

REDUCE is one of the older CASs (it was somewhat influenced by Macsyma). It is 
a commercial system. It is stable and efficient, and can solve larger problems in a 
given memory than, say, Mathematica. Its text-based interface looks old-fashioned; 
it has also an X interface, which provides 2-dimensional formula output (not of a 
very high quality) and a convenient on-line help. 

Here is a sample REDUCE session within T[T;Xmacs. 

Loading image file :/opt/reduce/lisp/psl/linux/red/reduce. img 
REDUCE 3.7, 15-Apr-1999, patched to 14-Jun-2001 ... 



1: (x"2-y"2)/(x~2+y~2)+sin(alpha)"2; 



sin (a) x^ -|- sin (a) y'^ 



X 



r 



2: (x+y-l)-5; 

x^ + bx^y - bx^ + 10x^2/^ - 20x^y + 10 x^ 
- 10 x^ 



.2„,3 



2 „,2 



lOx^r - 30x^2/^ + 30x^?/ 



bxy'^ — 20 x?/^ + 30 xt/^ — 20 xy + 5x + y^ 
- lOy"^ + 5y - 1 
3 : solve (a*x~2+b*x+c=0 , x) ; 
V - 4 a c + 62 _ 



X 



2a 



X 



(V-4ac + 62 + 6) 



2a 



4 : int (sqrt (x"2+a) , x) ; 

Va + x^x + log r ^^+^ 



5y* + 10 r 



5: int(exp(sin(x)) ,x) ; 

6: df(f(x),x,2); 
Declare f operator ? y 

d^ f jx) 

7: taylor(sin(x) ,x,0, 10) ; 
1 , 1 . 1 



X — - X 
6 



X 



120 5040 

8: m:=mat ((a,b) , (c,d)) ; 

fab 
9: 1/m; 

d -6 



x'^ + 



362880 



X 



O (x^i) 



ad — be ad — be 



ad — be ad — be 

10: plot(sin(x)/x,x=(-10 .. 10)); 

11: procedure f(x,y); sin(sqrt(x~2+y"2))/sqrt(x"2+y~2) ; 

/ 

12: plot(f (x,y) ,x=(-10 .. 10),y=(-10 .. 10) ,liidden3d,points=40) ; 

13: bye; 

Quitting 

The end 

Plots appear in separate windows, they are displayed by gnuplot (Fig. |]). Un- 
fortunately, it is not possible to control them interactively, one has to use options 
in the plot procedure. In order to save a plot to a PostScript file, the options 

terminal="postscript eps" ,output="f ilename.eps" 

are used. 



REDUCE Plot 





Figure 3: REDUCE plots 



The question mark icon displays help menu. REDUCE manual (written in I^T[t;X) 
is imported into TJ^macs. Some ETgX constructs are not handled correctly, but, 
nevertheless, the manual is quite readable. 

All interfaces described in this talk are preliminary, and require more work. It is 
not difficult to implement TEXmacs interfaces with more CASs. TgXmacs progresses 
fast; in the future, it can become a complete scientist's work place, suitable both 
for writing articles and for doing calculations using various external systems, within 
the same comfortable environment. 

I am grateful to Joris van der Hoeven for numerous discussions about TEXmacs 
and CAS interfaces; to William Schelter for his great help with Maxima and its KT^X 
generation; to Ralf Hillebrand for similar help with MuPAD and for providing an 
improved ETf^X generation library; to Winfried Neun for useful discussions about 
REDUCE. 
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